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Abstract 

Analytical representations in the time and frequency domains are de- 
rived for the most frequently used phenomenological fit functions for 
non-Debye relaxation processes. In the time domain the relaxation func- 
tions corresponding to the complex frequency dependent Cole-Cole, Cole- 
Davidson and Havriliak-Negami susceptibilities are also represented in 
terms of //-functions. In the frequency domain the complex frequency 
dependent susceptibility function corresponding to the time dependent 
stretched exponential relaxation function is given in terms of //-functions. 
The new representations are useful for fitting to experiment. 
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Analytical representaions of relaxation functions in the time domain and 
susceptibilities in the frequency domain arc important to fit experimental data 
in a broad variety of experiments on glasslike systems. Dielectric spectroscopy, 
viscoelastic modulus measurements, quasielastic light scattering, shear modulus 
and shear compliance as well as specific heat measurements all show strong 
deviations from the normalized exponential Debye relaxation function 

fit) = exp(-t/T) (1) 

where r is the relaxation time All relaxation functions in this paper are 
normalized to /(O) = 1. Relaxation in the frequency domain is described in 
terms of the normalized complex susceptibility 

^(u) = ^M^2^ {/(t)} (u) (2) 

Xo - Xoo 

where u = — 10, is the frequency, x(0) is a dynamic susceptibility normalized by 
the corresponding isothermal susceptibility, xo = liin0^o R-ex(0) is the static 
susceptibility, Xoo = lim0^oo Re x(0) gives the "instantaneous" response, and 
C{f{t)} (u) is the Laplace transform of the relaxation function f{t). For the 
exponential relaxation function this leads to 

1 + 10T 

i.e. the well known Debye susceptibility. 

Most generalizations of equations (0) and (||) for glasses and other complex 
materials are obtained by the method of introducing a fractional "stretching" 
exponent. In the time domain this method leads to the "stretched exponential" 
or Kohlrausch relaxation function given as 

f{t)^eMHt/r,f] (4) 

with fractional exponent /3 Q . Of course all formulae obtained by the method of 
stretching exponents are constructed such that they reduce to the exponential 
Debye expression when the stretching exponent becomes unity. Extending the 
method of stretching exponents to the frequency domain one obtains the Cole- 
Cole susceptibility Q 

the Davidson-Cole expression |^ 

or the combined Havrialiak-Negami form js] 

(1 + [uTHrP 

as empirical expressions for the experimentally observed broadened relaxation 
peaks. Most surprisingly, the analytical transformations between the time and 
frequency domain for general values of the parameters in these simple analytical 
expressions seem to be largely unknown m, and authors working in the time 
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domain usually employ the stretched exponential function while authors work- 
ing in the frequency domain use the stretched susceptibilities. An exception are 
the results in 0| where the real and imaginary part of the elastic modulus were 
obtained for Kohlrausch relaxation. Note however that there is a sign error 
in the real part in the results of |^]. It is therefore the purpose of this short 
communication to rederive expressions for the Kohlrausch susceptibility in the 
frequency domain. Secondly the same methods are used to obtain for the first 
time the relaxation function corresponding to the Havriliak-Negami suscepti- 
bility (and hence also the Cole-Davidson and Cole-Cole susceptibilities) in the 
time domain. It is hoped that these expressions will be useful for facilitating 
the fitting of experimental data. 

The objective of this paper is achieved by employing a method based on 
so called ff-functions The iJ-function of order {m,n,p,q) e and with 
parameters Ai e = 1, . . . ,p), Bi G R+(i = 1, . . . , g), G C(j = 1, . . . 

and 6i G C(i = 1, . . . , g) is defined for z G C, z 7^ by a contour integral in the 
complex plane m, M 



{ai,Ai), {ap,Ap) \ i 
(6i,i?i),...,(6„i3,) J2^^Jc 



77(s)z-" ds (8) 



where the integrand is 



Vis) = ^ . (9) 

Y[ r(a, + A,s) Yl r(l-fe, 

i—n~\-l 2— m+1 

In (|^) z"" = exp{— s log l^l — iargz} and argz is not necessarily the principal 
value. The integers m,n,p,q must satisfy 

< m < g, <n<p, (10) 

and empty products are interpreted as being unity. For the conditions on the 
other parameters and the path of integration the reader is referred to the lit- 
erature 1^ (see ||l^, p.l20ff] for a brief summary). The importance of these 
functions for glassy relaxation arises from the facts that (i) they contain most 
special functions of mathematical physics as special cases and (ii) their Laplace 
transform is again an i/-function. Moreover they possess series expansions that 
are generalizations of hypergeometric series. 

Based on the convenient properties of i7-functions it is possible to derive 
time and frequency domain expressions for all non-Debye relaxation functions 
and suceptibilities. The results of the calculations are summarized in the tables 
below. Table gives i7-function representations for all relaxation functions in 
the time domain. Table ^ gives the corresponding power series for all relaxation 
functions for small and large times. Table |^ summarizes i7- function represen- 
tations for the susceptibilities in the frequency domain, and Table ^ gives their 
power series expansions. Note that those power series where the domain of 
validity is given by a limit are asymptotic series. In these tables the notation 

r{a,x) = / y^'-'e-ydy (11) 
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denotes the complementary incomplete Gamma function, and the abbreviation 



oo ^. 



(12) 



is the Mittag-LefHer function. In addition the short hand notation 



11 
11 



( 



(1,1) 

(l,a) 



) 



(13) 



was introduced for writing the Kohh'ausch susceptibility. 

In summary the present paper has given unified representations of nonexpo- 
nential relaxation and non-Debye susceptibilities in terms of 77-functions. These 
representations lead to computable expressions that can be used to investigate 
the relations between the Kohlrausch susceptibility and other fit functions |jri| . 
It is hoped that the i/-function representations given here will help to facilitate 
the computational transformation between the frequency and time domain in 
theoretical considerations and experiment. 
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Table 1: iJ-function representation for the normalized relaxation functions 
(/(O) = 1) with relaxation time r. 
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Tabic 2: Scries expansions for normalized relaxation functions (/(O) = 1) with 
relaxation time r. Series are asymptotic whenever its range of validity is given 
as a limit. 
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Table 3: H- function representations for the normalized frequency dependent 
complex susceptibilities {u = —10). 
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Table 4: Scries representations for the normahzcd frequency dependent complex 
susceptibilities (it = —10). Series are asymptotic whenever its range of validity 
is given as a limit. 







series 


Debye 


1 

1 + UT 


00 

Y,{-lf{UTf \ut\ < 1 

fc=0 
00 

-Y^{-1)\UT)-^-' \UT\>1 


Kohlrausch 






Cole-Cole 


1 

1 + (UT„)« 


00 

^(-1)'=(UT„)"'= \UT^\ < 1 

k=0 
00 

-5^(-l)*=(iiT„)-«('=+i) |ur„|>l 

fc=0 


Cole-Davidson 


1 

(1 + UT^)-^ 


f(-i)'r(fc + 7). , 


Havriliak- Negami 


1 

(1 + {uthYY 





5 



